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ABSTRACT 


This  paper  presents  a nev  approach  to  the  analysis  of  finite  element 
methods  based  on  C® -finite  elements  for  the  approximate  solution  of  2nd  order 
bouniary  value  problems  in  which  error  estimates  are  derived  directly  in  terms 
of  two  mesh  dependent  norms  that  are  closely  related  to  the  norm  and  to 
the  2nd  order  Sobolev  norm,  respectively,  and  in  vdiich  there  is  no  assumption 
of  quasi -uniformity  on  the  mesh  family.  This  is  in  contrast  to  the  usual  analy- 
sis in  which  error  estimates  are  first  derived  in  the  1st  order  Sobolev  norm  and 
subsequently  are  derived  in  the  L^  norm  and  in  the  2nd  order  Sobolev  norm  - the 
2nd  order  Sobolev  norm  estimates  being  obtained  under  the  assumption  that  the 
functions  in  the  underlying  approximating  subspaces  lie  in  the  2nd  order  Sobolev 
space  and  that  the  mesh  family  is  quasi-uniform. 
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SIGNIFICANCE  AND  EXPLANATION 

Over  the  last  twenty  years  there  has  been  an  extensive  development  of  finite 
element  methods  for  the  approximate  computer  solution  of  problems  in  continuum 
mechanics  (fracture,  combustion,  fluid  flow,  etc.).  The  exact  solution  u that 
we  are  looking  for  satisfies  a differential  equation  boundary  value  problem.  We 
seek  an  approximation  u^  to  u in  a well-defined  family  of  approximating  func- 
tions S . For  example,  S might  consist  of  continuous,  piecewise  linear  func- 
n n 

tions,  there  being  n subregions  in  each  of  which  any  member  of  is  linear. 

• i 

We  would  expect  that  the  larger  n is,  the  closer  u^  is  to  u. 

The  key  type  of  result  used  in  the  theoretical  justification  of  such  methods 

is  an  estimate  of  the  form 

(•)  ||u  - u^ll  < c Inf  ||u  - xll  • 

Here  the  norm  ||u  - vj|  denotes  some  measure  of  the  difference  between  u and 
V.  The  right  hand  side  is  a constant  (independent  of  n)  times  the  lower  bound 
of  the  difference  between  u and  any  X ii^  This  estimate  reduces  the  pro- 

blem of  determining  the  error  in  u^  to  a much  easier  one  in  approximation  theory; 
namely,  how  closely  can  u be  approximated  by  a X in  S^. 

There  are  several  non-ideal  features  of  the  existing  theory  — for  certain 
norms,  estimates  of  the  type  (*)  are  not  known  (in  many  situations  estimates  of 
the  form  ||u  - u^|l  j<  C h'*  are  known,  where  h is  the  size  of  the  largest  mesh 
subdomain  and  a is  some  positive  constant) , and  for  certain  norms  quasi-uniform 
meshes  are  required  (meshes  for  which  the  ratio  of  maximum  to  minimum  mesh  size  is 
bounded  independent  of  the  maximum  mesh  size) . As  a consequence,  local  mesh  re- 
finement is  not  accounted  for.  The  present  paper  presents  an  approach  to  such 
methods  in  which  (a)  estimates  of  the  type  (*)  are  obtained,  (b)  the  results  are 
valid  for  arbitrary  meshes,  and  (c)  local  mesh  ref iuc-ment  is  properly  accounted  for. 

The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive  summary 
lies  with  MRC,  and  not  with  the  authors  of  this  report. 
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Then  let  v (x)  » j f(t)dt.  v is  clearly  in  S..  Now,  using  the  Braable*>Kilbert 
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In  Section  6 we  apply  the  results  of  Section  2 with 
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